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MECHANICS  OF  COMPOSITE  MATERIALS 


1.  RESEARCH  OBJECTIVES 

Research  under  the  subject  Contract  represents  a continuing  effort3 
to  construct  viable  nonlinear  micro-mechanical  continuum  models  to  pre- 
dict the  thermomechanical  response  of  advanced  composite  materials.  The 
thermomechanical  processes  under  study  include: 

. Linear  and  nonlinear  wave  propagation 

. Linear  and  nonlinear  diffusion  (thermal  and  hygrothe rmal) 

. Mechanisms  and  accumulation  of  damage. 

The  thermomechanical  measures  include: 

. Global  stress,  deformation,  temperature,  internal  energy  fields 
. Local  stress,  deformation,  temperature,  internal  energy  fields 
. Residual  properties  of  damaged  composites. 

2.  RESEARCH  PROGRESS 

2.  1 Discussion 

Excellent  progress  has  been  made  to  date  in  modeling  both  laminated 
and  fibrous  composites  as  homogeneous  continua  with  microstructure.  The 
resulting  theory  is  generally  referred  to  as  a "Theory  of  Interacting  Continua" 
or  a "Mixture  Theory  with  Microstructure".  In  contrast  to  classical  theories 

3Previous  work  was  conducted  under  AF-AFOSR  Grant  70-1975  and  AF-AFOSR 
Contract  75-2870. 


knowledge  of  the  initial  microstructural  geometry  of  the  composite,  the 
component  constitutive  relations,  and  the  component  interface  physics. 


In  further  contrast  to  classical  theories,  the  current  model  provides  infor- 
mation on  stress,  displacement,  strain,  temperature,  and  internal  energy 
fields  within  the  micro  components  as  well  as  global  measures  of  these 
quantities.  The  present  theory  also  accounts  for  debonding  at  the  micro- 
component interfaces  - a basic  damage  mode.  Models  incorporating  crack 
distributions  throughout  the  constituents  are  under  study. 

The  most  general  mathematical  model  considered  to  date  is  a non- 
lineai , anisotropic,  thermomechanical  mixture  theory  with  microstructure. 
Particular,  simplified  farms  of  this  model  have  been  developed  to  cover 
special  cases  of  composite  geometry,  material  constitutive  relations,  and 
physical  processes  or  problems.  In  all  cases,  however,  theoretical 
construction  is  based  upon  one  of  two  basic  methods  which  are  in  turn 
based  upon  an  asymptotic  scheme  in  which  dominant  signal  wavelengths 
associated  with  the  physical  process  are  assumed  to  be  large  in  comparison 
with  typical  composite  microdimensions.  An  explanation  of  the  foregoing 
construction  methods,  as  well  as  a number  of  interesting  examples,  can 
be  found  in  Appendix  II. 

The  composites  studied  using  the  above  asymptotic  technique  range 
from  elementary  laminated  to  complex  three-dimensional  advanced 
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composites.  In  addition  to  certain  model  composites,  materials  have 
included  quartz  phenolic,  carbon  phenolic,  carbon  carbon,  and  graphite 
epoxy.  Representative  problem  areas  include  nose  tips,  heat  shields, 
engine  blades,  nozzles,  and  hygrothermal  environments. 

Cases  completed  and  in  publication  thus  far  include: 

. A general  theory  for  wave  propagation  in  linear  elastic  and 
linear- viscoelastic  laminated  composites  with  isotropic  laminae; 

. A general  theory  for  linear  thermal  diffusion  in  laminated 
composites  with  isotropic  laminae; 

. A "first  order"  theory  for  linear  wave  propagation  in  anisotropic 
laminates  with  oblique  layup; 

. "First  order"  theories  for  wave  propagation,  diffusion,  and 
debonding  in  unidirectional  fibrous  composites  with  cylindrical  fibers  of 
arbitrary  cross  section; 

. "First  order"  nonlinear,  anisotropic  and  isotropic,  thermo- 
dynamic theories  for  wave  propagation,  diffusion,  and  damage  accumulation 
due  to  debonding  in  both  laminated  and  fibrous  composites  of  specific  types. 

In  the  general  linear  cases  above,  a hierarchy  of  models  is  system- 
atically defined  by  the  order  of  the  truncation  of  the  asymptotic  sequence 
involved  in  the  theoretical  construction  process.  Retention  of  all  terms 
in  this  sequence  leads  to  a formally  exact  theory.  Retention  of  only  the 
zeroth-order  plus  first  order  terms  results  in  a so-called  "first  order" 


theory  which  can  be  cast  in  the  form  of  a mixture  theory.  The  latter 
usually  constitutes  an  adequate  micromechanical  model  of  a composite 
for  engineering  purposes.  All  nonlinear  models  are  "first  order"  theories. 

The  first  order  models  for  laminates  and  fibrous  composites 
incorporate  the  effects  of  arbitrary  fiber  geometry,  anisotropy,  and  inter- 
face debonding.  In  addition  to  diffusion,  a thermodynamic  theory  of  elastic- 
plastic  wave  propagation  under  finite  deformation  has  been  developed.  The 
constitutives  laws  of  the  microcomponents  currently  in  publication  relate 
the  mean  stress  (trace  of  the  stress  tensor)  to  density  and  internal  energy 
by  a Mie-Griineisen  caloric  equation  of  state  and  the  stress  deviator  ten- 
sor to  <-he  rate-of-deformation  tensor  by  an  elastic-perfectly  plastic 
relation  of  the  von  Mises  type;  extension  to  general  anisotropic  plasticity 
utilizing  a plastic  potential  and  an  associated  flow  ride  is  under  develop- 
ment. Constituent  debonding  and  component- interface  physics  is  governed 
by  a Coulomb  frictional- type  relation;  the  latter  constitutes  the  present 
measure  of  "damage";  extension  of  this  measure  to  include  microcracking 
throughout  the  components  is,  as  was  noted  previously,  under  study. 

In  the  case  of  linear  elastic  wave  propagation,  model  accuracy 
has  been  demonstrated  by  comparison  of  phase  velocity  spectra  with  exact 
and/or  experimental  data.  Accuracy  superior  to  existing  theories  has 
been  observed,  as  well  as  first  order  model  applicability  down  to  wave- 
lengths on  the  order  of  typical  composite  microdimensions  (when  the  energy 

« 

I 

4 

— i 


is  partioned  primarily  in  the  first  mode).  In  addition,  a number  of  transient 
pulse  cases  have  been  treated;  these  exhibit  good  correlation  with  exact  and/or 
experimental  data.  A sampling  of  typical  transient  pulse  results  for  both 
laminated  and  unidirectional  fibrous  composites  is  illustrated  in  Figs.  1-4 
of  Appendix  1.  Typical  transient  wave  propagation  results  for  a complex, 
advanced,  three-dimensional  composite  are  given  in  Figs.  8-10  of  Appendix  I. 
Complete  details,  including  extensive  phase-velocity-spectra-calculations, 
can  be  found  in  the  appropriate  references  listed  under  Section  2.  3. 

In  the  case  of  large  deformations  and  elastic-plastic  wave  propaga- 
tion, theoretical  transient  pulse-results  for  propagation  parallel  to  the 
fibers  of  a unidirectional  fibrous  composite  have  been  compared  with  data 
from  a 2D-Lagrangian  finite  difference-code  (CRAM);  the  latter  is  essentially 
exact.  As  can  be  observed  from  Appendix  I,  Figs.  11-14,  the  agreement 
between  "exact"  and  approximate  calculations  was  found  to  be  excellent. 

Model  accuracy  in  the  case  of  interface  debonding  has,  to-date, 
been  assessed  by  comparison  with  experimental  data  on  delaminated 
plates  composed  of  alternating  layers  of  Polymethyl  Methacrylate  (PMMA, 
Rohm  and  Haas  Type  A),  and  6061-T6  aluminum.  The  laminae  of  the 
composite  were  oriented  perpendicular  to  the  impact  plane,  and  struck 
by  a projectile  fired  from  a 10  cm  bore  light  gas  gun.  On  the  basis  of 
averaged  particle  velocity  measured  at  the  rear  face  of  the  specimen, 
the  model  hypothesized  yields  excellent  correlation  with  the  experimental 
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results.  The  test  setup  and  typical  results  are  illustrated  in  Figs.  15-18 
of  Appendix  I. 

In  addition  to  the  above,  studies  concerning  the  accuracy  of  the 
thermal-portion  of  the  continuum  models  were  undertaken.  In  particular, 
diffusion-type  processes  in  a laminated  composite  with  periodic  micro- 
structure were  examined.  Solutions  for  the  lowest-order  models  were 
compared  with  "exact"  results  from  a finite  difference  code.  For  most 
cases  the  lowest-order  "effective  conductivity"  theory  was  found  to  yield 
good  results.  For  exceptional  problems  requiring  a higher-order  theory, 
a modified  version  of  the  lowest-order  theory  was  found  to  yield  excellent 
correlation  between  exact  and  approximate  solutions.  For  man.  problems 
of  the  diffusion- type,  these  elementary  equations  may  offer  an  attractive 
alternative  to  other  means  for  obtaining  solutions.  Typical  results  for 
the  case  of  heat  propagating  normal  to  the  layers  of  the  laminated  com- 
posite of  Fig.  1,  Appendix  I,  are  illustrated  in  Figs.  19  and  20  of 
Appendix  I.  In  each  case  a uniform  temperature  Tq  is  applied  to  y = 0 
from  t = 0 and  t = tQ  while  the  boundary  temperature  is  zero  for  all  other 

time.  (Here  A = h,  + h„,  p c = heat  capacity  of  a-material,  k = 

1 2 ra  or  a 

thermal  diffusivity).  It  is  emphasized  here  that,  although  Figs.  19,  20 
refer  to  a one-dimensional  example,  the  theory  is  valid  for  three- 
dimensional  problems.  The  one-dimensional  case  was  selected  mainly 
because  the  most  severe  test  of  the  theory  occurs  when  heat  flow  is 
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normal  to  the  laminates  (since  material  properties  are  discontinuous  in 
this  direction). 

In  the  last  interim  scientific  report,  it  was  noted  that  a major  develop- 
ment in  theoretical  modeling  was  effected  during  the  previous  research  period. 
The  development  allows  treatment  of  arbitrary  fiber  geometry  and  arrays. 
Basically,  the  task  of  constructing  mixture  theories  has  been  partioned  into 
two  isolated  steps.  The  first  step  involves  the  determination  of  the  general 
mixture  equations  for  a particular  problem  class;  these  equations  involve  coef- 
ficients which  must  be  determined  from  the  detailed  microstructural  geometry, 
interface  data,  and  constitutive  relations.  The  second  step  defines  the  constants 
via  a quasi-static  "microboundary  -alue  problem"  (MVP). 

Resolution  of  the  MVP  is  carried  out  by  a variational-based  finite  element 
method.  The  procedure  has  been  applied  to  diffusion  and  wave  propaga- 
tion in  unidirectional  composites.  Here  important  information  has  been 
obtained  regarding  the  accuracy  of  the  concentric  cylinders  approximation 
used  frequently  in  practice  for  circular  fibers  in  a hexagonal  array  and 
for  rectangular  fibers  in  a similar  array.  Extension  of  the  procedure  to 
arbitrary  three  dimensional  geometries  is  under  development. 

A discussion  of  the  above  subject  matter  can  be  found  in  the  survey 
paper  listed  under  Section  2.  3.  For  the  readers  convenience,  the  latter 
is  included  as  Appendix  II  of  this  report. 
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Finally,  during  the  current  research  period  another  major  advance 


in  our  theoretical  modeling  procedure  was  accomplished.  This  development 
allows  a unified  treatment  of  laminated  and  fibrous  composites,  and  provides 

a substantial  increase  in  model  accuracy.  The  mixture  theory  construction 

" 

procedure  in  this  case  is  based  upon  the  Regular  Asymptotic  Method  in  conjunc- 
tion with  the  use  of  a multi- variable  expansion  procedure.  The  approach  has 
been  successfully  applied  to  both  wave  propagation  and  heat  conduction  in  lami- 
nated composites.  Fibrous  composites  are  under  study.  The  resulting  form 
of  the  mixture  theory  indicates  a significant  flaw  currently  exists  in  some 
general  phenomenological  mixture  theories.  Reports  and  publications  on  this 
subject  are  in  preparation. 

2.  2 Research  Applications 

The  fundamental  concepts  embodied  in  the  aforementioned  wave 
propagation  studies  have  been  successfully  employed  in  an  applied  research 
effort  to  model  1-D  (FMI  COMRAD)  and  AVCO  3DQP  quartz  pheon’ic 
composites.  These  models,  and  the  resulting  numerical  code  (TINC)  are 
presently  in  practical  use  by  the  engineering  profession.  Typical  results 
for  3 DQP  a-e  illustrated  in  Figs.  5 - 10  of  Appendix  I. 

2.  3 List  of  Publications 

The  following  publications  are  representative  of  the  research 
progress  attained  during  previous  and  present  Grant  periods.  All 
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APPENDIX  I 


Compendium  of  Typical  Theoretical  Results 

Presented  in  this  Appendix  is  a sampling  of  typical  transient  pulse 
and  diffusion  calculations  carried  out  using  the  Theory  of  Interacting 
Continua  (TINC).  The  cases  discussed  include  laminated  and  unidirectional 
fibrous  composites  as  well  as  a complex  3-dimensional  layup.  Results  of 
the  TINC  theory  for  each  example  are  compared  to  experimental  data. 


Laminated  Composites 

The  geometry  of  the  laminated  composites  studies  is  illustrated 
in  Fig.  1.  One  problem  of  interest  is  that  of  symmetric  (P)  wave 
propagation  parallel  to  the  layers,  with  material  behavior  restricted  to 
the  linear  elastic  range.  A comparison  of  TINC  and  experimental  results 
is  given  in  Fig.  2.  The  experimental  results  were  reported  by  Whittier 
and  Peck  [1]  for  specimens  composed  of  Thornel  (high  modulus  carbon) 
fibers  reinforcing  a carbon  phenolic  matrix.  Specimens  of  1/4-inch 
thickness  were  subjected  to  a uniform  pressure  at  the  left  boundary,  with 
a step  function  time-dependence  induced  by  a gas  dynamic  shock  wave  of 
about  70  psi.  The  numerical  calculations  were  initiated  by  impacting  a 
step  function  velocity  of  7.786  cm/sec  to  both  constituents  at  the  boundary. 
While  this  condition  does  not  correspond  precisely  to  the  experiment,  it 
was  felt  that  the  error  introduced  would  be  negligible  when  far  removed 
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from  the  front  surface.  Figure  2 depicts  the  comparison  between 
experiment  and  theoretical  results  for  averaged  rear  surface  velocity, 
normalized  on  the  boundary  input  velocity.  Since  absolute  times  were 
not  measured  in  the  tests,  theoretical  and  experimental  results  were 
matched  at  their  respective  first  peak  arrivals.  In  the  figure,  results 
from  the  full  TINC  equations  solved  numerically  are  denoted  "Continuum 
mixture  theory.  " The  results  denoted  "Simplified  theory"  are  for  an 
elementary  form  of  the  TINC  theory  admitting  a closed-form  solution. 

Unidirectional  Fibrous  Composite 

The  TINC  theory  has  been  applied  to  the  problem  of  wave  prop- 
agation parallel  to  the  fibers  of  a unidirectionally  reinforced  composite 
with  hexagonal  array.  The  hexagonal  array  is  approximated  by  concentric, 
linearly  elastic  cylinders  as  illustrated  in  Fig.  3. 

In  Fig.  4,  transient  pulse  data  predicted  by  TINC  is  compared 
with  experimental  data  on  a unidirectional  quartz  phenolic  fibrous 
composite.  The  experiments  were  subjected  to  a 70  psi  step  function 
in  pressure  via  a shock  tube.  The  figure  depicts  the  comparison  between 
the  experimental  and  TINC  theory  code  predictions  of  area-averaged  rear 
surface  velocity:.  Results  of  the  simplified  TINC  theory  are  also  por- 
trayed. Once  again  the  experimental  and  theoretical  results  are  matched 
at  their  first  peak  arrivals  since  absolute  arrival  times  were  not  measured 
experimentally.  As  was  the  case  in  Fig.  2,  the  times  shown  are  those 
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predicted  by  the  theory. 

Three-dimensional  Layups 

Predictions  of  transient  pulse  behavior  from  the  TINC  theory 
have  been  compared  to  experimental  results  for  a 3-dimensional  quartz 
phenolic  composite  (3DQP)  manufactured  by  AVCO.  A representative 
selection  of  the  results  are  given  here. 

A "block"  3DQP  specimen  (Fig.  5)  was  subjected  to  shock  tube 

tests  by  the  Aerospace  Corporation.  The  test  procedure  is  the  same 

as  that  used  for  laminated  and  fibrous  composites  and  is  described  above. 

Once  again  the  theoretical  and  experimental  boundary  conditions  are  not 

identical,  but  this  decrepancy  is  expected  to  introduce  little  error.  The 

averaged  rear  surface  velocity  measured  by  experiment  is  compared  to 

that  predicted  by  the  TINC  theory  in  Fig.  6.  Also  shown  are  results 

predicted  by  the  elementary  "Head  of  the  Pulse"  approximate  [2]. 

The  experimental  data  has  been  forced  to  coincide  with  the  theoretical  at 

V/V  =1/3  since  absolute  arrival  times  were  not  measured, 
o 

Another  set  of  relatively  low  pressure  wave  propagation  experiments 
were  performed  at  the  Air  Force  Weapons  Laboratory  using  a light  gas 
gun  and  "midspace"  3DQP  specimens,  shown  in  Fig.  7.  In  this  test,  three 
specimens  of  different  thicknesses  were  mounted  on  a single  target  and 
impacted  with  a .01-inch  mylar  flyer.  The  input  stress  profiles  for  the 
three  specimens  are  illustrated  in  Fig.  8.  For  the  TINC  calculations. 
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input  velocity  time-histories  corresponding  to  the  stress-time  data  of 
Fig.  8 were  used.  A comparison  of  experimental  and  theoretical  results 
for  the  thickest  of  these  specimens  is  shown  in  Fig.  9.  Since  absolute 
arrival  times  were  not  measured,  the  two  curves  are  forced  to  coincide 
at  their  peak  pressures.  The  agreement  here  is  excellent,  especially 
in  view  of  the  fact  that  the  wave  is  measured  only  2-1/2  microdimensions 
from  the  front  surface.  Also,  with  the  10kb  input  stress,  this  is  a non- 
linear calculation. 

The  AFWL  also  performed  a series  of  high-pressure  gas  gun 
experiments  on  "midspace"  3DQP  specimens  (Fig.  7).  Experimental 
and  theoretical  data  are  compared  r«  Fig.  10.  The  observed  agreement 
indicates  that  the  thermodynamic  and  nonlinear  constitutive  behavior 
incorporated  into  the  TINC  theoretical  model  are  adequate.  Note  that 
here  the  rear  surface  velocity  is  measured  only  1-1/2  microdimensions 
from  the  front  surface.  Thus,  as  before,  the  TINC  theory,  though 
derived  as  an  "outer"  solution,  yields  surprisingly  good  results  even 
at  locations  quite  near  the  boundary. 

Comparison  of  TINC  and  2-D  Finite  Difference  Code 

It  is  of  interest  to  compare  the  results  from  a TINC  calculation 
to  those  of  a typical  two-dimensional  Lagrangian  finite  difference  code 
which  is  essentially  exact.  Such  a code  (CRAM)  is  available  at  Systems, 

Science  and  Software,  Inc.  The  two  codes  were  compared  for  the 
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concentric  cylinder  geometry  of  Fig.  11.  Computer  costs  for  CRAM  were 
approximately  500  times  those  for  running  the  TINC  code.  Several  typical 
results  are  illustrated  in  Figs.  12-14  for  a square  wave  input  velocity 
(shown  as  problem  #2  on  Fig.  11).  Once  again,  the  calculations  were 
carried  out  within  a few  microdimensions  of  the  boundary. 

Debonding  of  Constituent  Interfaces 

Model  accuracy  in  the  case  of  interface  debonding  or  cracking 
has,  to  date,  been  assessed  by  comparison  with  experimental  data 
dealing  with  delaminated  plates  composed  of  alternating  layers  of 
Polymethyl  Methacrylate  (PMMA,  Rohm  and  Haas  Type  A),  and  6061-T6 
aluminum.  The  laminae  of  the  composite  were  oriented  perpendicular 
to  the  impact  plane,  and  struck  by  a projectile  fired  from  a 10  cm  bore 
light  gas  gun.  On  the  basis  of  averaged  particle  velocity  measured  at 
the  rear  face  of  the  specimen,  the  model  hypothesized  yields  excellent 
correlation  with  the  experimental  results. 

The  test  setup  employed  is  illustrated  in  Fig.  15.  Typical 
correlation  between  theoretical  and  experimental  results  is  shown  in 
Figs.  16-18.  Complete  details  can  be  found  in  reference  [13^ 

Diffusion 

The  laminated  geometry  of  Fig.  1 was  selected  to  evaluate  the 
continuum  model  construction  procedure  for  diffusion  problems  (a 
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complete  discussion  of  this  work  can  be  found  in  reference  [4]).  The 
case  of  heat  propagation  normal  (in  the  y-direction)  to  the  layers  was 
selected  for  a comparison  of  exact  and  continuum  theory  solutions  since 
it  is  known  that  this  case  represents  the  most  severe  test  of  the  latter. 
Figure  19  represents  a case  of  similar  material  properties  whereas 
Fig.  20  illustrates  the  comparison  for  widely  differing  material  prop- 
erties. Continuum  theory- results  were  obtained  via  closed-form 
solutions  for  a square-wave  temperature  input  at  the  boundary  y = 0. 
"Exact"  results  were  obtained  using  a finite  difference  code.  The 
calculations  illustrate  the  ability  of  the  first-order  theory  to  model 
micrc-xtructure  details  and  closely  match  the  exact  data. 
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Fig.  2.  Comparison  of  TINC  and  experimental  rear  surface  velocity 
results  for  thornel  reinforced  carbon  phenolic  laminates. 


Comparison  of  TINC  and  experimental  rear  surface 
velocity  results  from  quartz  fiber  reinforced  phenolic 
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Fig.  9.  Stress  history  at  rear  of  0.  254  in.  3DQP 
plus  0.060  in.  plexiglass  (shot  326). 
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Fig.  11.  Geometry  for  2D  CRAM  Calculation. 
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Fig.  12.  Comparison  of  TINC  and  average  CRAM  particle 
velocity  histories  for  radial  fiber. 
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Fig.  13. 
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Comparison  of  TINC  and  average  CRAM  particle 
velocity  histories  for  sheath. 
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Fig.  14.  Comparison  of  TINC  and  CRAM  predictions  of 
axial  stress  (x  = 0.  2 cm,  t = 1.0  psec). 
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Fig,  15.  Experimental  arrangement  for  Sandia  experiments. 
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Fig.  16.  Comparison  of  theory  with  experiment  No.  2. 
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Fig.  17.  Comparison  of  theory  with  experiment  No.  6. 
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Fig.  20.  Temperature  distribution  at  r = 1 for  widely 

differing  material  properties  and  small  t . 
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mixture  theories  with  microstructure  for  wave 

PROPAGATION  AND  DIFFUSION  IN  COMPOSITE  MATERIALS 

Two  procedures  for  constructing  mixture  theories  with  micro- 
structure are  discussed  for  wave  propagation  and  thermal  diffusion 
processes  in  binary  laminated  and  fibrous  composites.  Each  method 
proceeds  from  the  composite  microstructure  and  is  based  upon  an 
asymptotic  scheme  wherein  the  ratio  of  trans verse-to-longitudinal 
characteristic  times  associated  with  a physical  process  is  assumed  to  be 
small.  The  methods  retain  considerable  information  regarding  the 
mechanical  and/or  diffusive  fields  within  the  microcomponents. 

Theoretical  results  are  compared  with  both  experimental  and  exact  data 
in  an  effort  to  evaluate  model  accuracy.  The  presentation  includes  treat- 
ment of  geometric,  constitutive,  and  interface  nonlinearities. 

1.  INTRODUCTION 

A considerable  number  of  continuum  theories  have  been  proposed 
for  thermomechanical  processes  in  composites.  The  continuum  theory 
of  mixtures  represents  one  of  the  more  important  and  successful  theoreti- 
cal descriptions  of  such  multiphase  materials.  According  to  this  concept, 
the  composite  constituents  are  modeled,  at  each  instant  of  time,  as  super- 
posed continua  in  space.  Each  continuum  is  allowed  to  undergo  individ- 
ual deformations.  The  microstructure  of  an  actual  composite  is  then 
simulated  by  specifying  the  interactions  between  the  continua. 

Mixture  theories  for  multiphase  materials  have  been  developed 
from  both  macrostructural  and  microstructural  viewpoints.  The  vast 
majority  of  activity,  however,  has  been  directed  toward  the  former, 
which  concerns  the  formulation  of  a general  theoretical  framework  for 
classes  of  mixtures  based  upon  certain  phenomenological  postulates  that 
avoid  detailed  microstructural  considerations.  Excellent  accounts  of 
this  approach  can  be  found  in  [1-3]. 

The  development  of  mixture  theories  from  the  microstructural 
viewpoint  represents  a more  difficult  task.  The  reward,  however,  is  a 
model  - which  shall  be  called  a mixture  theory  with  microstructure  - 
that  is  completely  determined  by  the  geometry  and  constitutive  relations 
of  the  individual  components  and  which  provides,  to  a certain  degree  of 
accuracy,  information  concerning  microstruct ure  fields.  Of  course, 
the  price  one  must  pay  for  such  detail  is  often  a case-by-case  study 
rather  than  a general  theory.  Nevertheless,  it  should  be  recognized 
that,  while  general  macrostructural  - based  theories  may  appear  more 
elegant,  such  theories  may  not  be  practical  for  many  applications  due  to 
the  unreasonable  burden  placed  upon  the  experimentalist  to  determine  a 
multitude  of  constants  and/or  functions. 
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In  this  presentation,  focus  is  placed  upon  several  recent  works  by 
the  author  and  co-workers  concerning  the  construction  of  mixture  theories 
with  microstructure  for  laminated  and  fibrous  composites.  The  discus- 
sion is  limited  to  thermomechanical  processes  in  composites  with  initial 
periodic  microstructure.  Particular  emphasis  is  placed  upon  wave 
propagation  and  thermal  diffusion. 

The  presentation  is  divided  into  five  sections.  The  methodology 
associated  with  two  theoretical  construction  techniques  is  given  in 
Section  2 along  with  a brief  survey  of  applications.  These  techniques 
are  demonstrated  via  two  examples  in  Sections  3 and  4.  Finally,  in  an 
effort  to  illusidate  the  applicability  and  utility  of  the  techniques,  two 
applications  to  nonlinear  problems  are  discussed  in  some  detail  in 
Section  5;  these  include  comparisons  between  mixture  theory,  exact  and/ 
or  experimental  results. 

2.  METHODOLOGY  AND  SURVEY 

A variety  of  construction  procedures  may  be  used  to  generate 
mixture  theories  with  microstructure.  Two  such  procedures  are  selected 
for  discussion  herein  in  view  of  their  observed  success  in  describing 
thermomechanical  processes  in  composite  materials. 

The  first  procedure,  which  is  illustrated  in  Section  3 by  an 
elementary  example,  is  the  Differential-difference  Method.  This 
approach  has  been  used  to  develop  continuum  models  of  wave  propagation 
in  linear  elastic  [4-9,26]  and  linear  viscoelastic  [9-10]  laminated  com- 
posites, wave  propagation  in  linear  elastic  fibrous  composites  [7-11], 
wave  propagation  and  debonding  in  laminated  composites  [12],  and 
thermal  diffusion  in  laminated  composites  [13,  14],  The  method  consists 
of  eight  basic  construction  steps  which  may  be  summarized  as  follows: 

1)  the  field  quantities  (all  dependent  variables)  of  each  composite 
constituent  are  expanded  in  a suitable  spatial  series  (this  expansion, 
which  amounts  to  discretization  of  the  component  domains,  may  range 
from  a Taylor  series  for  sufficiently  simple  geometries  to  topological 
multiplexes  [15],  e.  g.  , finite  elements,  for  complex  geometries); 

2)  recurrence  relations  for  the  expansion  coefficient-functions  are  ob- 
tained by  use  of  the  field  equations  and  the  series  is  telescoped;  3)  applica- 
tion of  interface  continuity  conditions  lead  to  differential  (in  time)  - 
finite  difference  (in  space)  equations;  4)  an  associated  set  of  differential 
(in  time)  - functional  difference  equations  is  observed  to  contain,  as  a 
solution  subset,  solutions  of  the  finite  difference  relations;  5)  solutions 

of  the  functional  equations  are  smooth  fields  which  take  on  exact  values 
at  discrete  spatial  points;  these  smooth  fields  are  now  employed  as 
dependent  variables;  6)  all  spatial  differences  are  expanded  in  a Taylor 
series  - an  operation  which  furnishes  partial  differential  equations,  i.  e.  , 
a continuum  model;  7)  nondimensional  variables  and  a small  parameter 
representing  the  ratio  of  typical  micro-to-macro  dimensions  is  introduced 
and  the  operators  are  truncated  according  to  an  asymptotic  procedure; 

8)  mixture  equations  result  by  algebraic  manipulation. 

The  second  procedure,  illustrated  in  Section  4 by  an  elementary 
example,  is  the  Regular  Asymptotic  Method.  This  approach  has  been 
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used  to  model  nonlinear  wave  propagation  in  laminated  [16]  and  fibrous 
[17-19]  composites,  and  thermal  diffusion  in  fibrous  composites  [20-21], 
The  method  consists  of  six  basic  construction  steps  which  may  be  sum- 
marized as  follows:  1)  based  upon  the  particular  composite  and  the 
process  considered,  an  estimate  of  the  order  of  magnitude  of  all  field 
variables  is  made;  2)  the  field  equations  are  re-scaled  such  that  all 
dependent  variables  are  0(1);  in  the  process  a small  parameter  is  intro- 
duced which  constitutes  a measure  of  the  ratio  of  micro-to-macro- 
dimensions  of  the  problem;  3)  the  conservation  equations  are  averaged 
to  obtain  a standard  mixture  form;  4)  the  scaled  field  variables  are 
expanded  in  a regular  asymptotic  series;  5)  the  two  lowest  order  systems 
(including  interface  conditions)  are  used  to  obtain  expressions  for  the 
interaction  terms  and  the  constitutive  relations;  6)  the  constants  involved 
in  5)  are  obtained  from  the  solution  of  a static  micro  boundary  value 
problem  defined  over  a unit  cell;  a variational  principle  - based  finite 
element  method  is  suggested  to  resolve  this  boundary  value  problem. 


A judicious  choice  of  a particular  method  depends  upon  the 
problem.  The  Differential-difference  Method  is  best  suited  to  linear 
problems  and  elementary  geometries  such  as  laminates.  Here  formally 
exact  relations  may  be  obtained  in  many  cases.  The  Regular  Asymptotic 
Method  is,  on  the  other  hand,  applicable  to  nonlinear  problems  and 
complex  geometries.  In  this  case  one  is  usually  content  with  the  first 
few  terms  of  the  asymptotic  series. 


3.  EXAMPLE  1:  WAVE  PROPAGATION 


For  the  purpose  of  illustrating  the  Differential-difference  Method, 
consider  wave  propagation  in  a periodic  array  of  linearly  elastic,  homo- 
geneous and  isotropic  laminae,  perfectly  bonded  at  all  interfaces  as 
illustrated  in  Fig.  1.  Let  y = X£  for  notational  convenience  and  assume 
a state  of  uniaxial  motion  in  the  y-direction;  i.  e.  , the  case  of  one- 
dimensional longitudinal  waves  propagating  normal  to  the  laminae.  Re- 
stricting attention  to  small  isothermal  deformations,  the  appropriate 
conservation  and  constitutive  equations  are 


. 2 .(tt,k) 

OyO  - p&t  u) 


(cr-(x+2p)aju)(0'’k)  = 0 


(1) 


while  the  interface  continuity  conditions  are 

^ r nl (2,k)  (2)  . 

[u,o]  (h  ,t)  = [u,0]  (-h  ,t)  , 


wP-VU 


(2) 


In  the  above  the  superscripts  (tt,k)  refer  to  the  k layer  of  the  Q 
constituent,  a = 1 , 2;  k = 1 , 2,  3.  . . . The  variable  y'^A)  j.s  a local 

coordinate  with  origin  at  the  midplane  of  the  (a,k)  layer;  0 h 022  anc* 
u = U£are  normal  stress  and  displacement  (subscripts  on  ct  and  u have 
been  dropped  for  brevity).  The  quantities  p(cO,  h(<x),  \(°0 » are  mass 

density,  layer  half  thickness  and  Lame'  constants  for  the  a - constituent. 
In  addition  >y(  ) =■  g(  )/?y,  8t(  ) = 8(  )/gt. 

Wave  motion  in  the  composite  is  completely  specified  by  (1),  (2) 
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together  with  initial  conditions  at  t=0  and  boundary  data  on  y=0,  H(H  s «). 
The  objective  of  the  subsequent  work  is  to  replace  this  set  of  equations 
(infinite  if  H = =)  by  a single  set  of  partial  differential  equations  which 
represent  a continuum  mixture  theory. 

3. 1 Expansion,  Recurrence  Relations 

Let  us  expand  the  stress  and  displacement  about  each  lamina  mid- 
plane as  follows: 

,n 

(3) 


, (Or,k)  - ...  (Or.k)'  , 

g(y  ,t)  = I g,  ,(t)y  /ni 

n=0  (n) 


(ot  k)  (q 

Here  g represents  any  of  o , u . Upon  substituting  (3)  into  (1)  and 
equating  powers  of  y^a>k),  one  obtains  the  differential  - recurrence 
relations 

.(«.k)  _ 2,  ,(o,k)  (a,k)  -r(al  (a,k) 

C ■ %>  =E  V+l)  • <4a) 


Here 


(a)  ~(a).  (a)  ~(a)  (a) 

c = E /p  , E = (X+2  (a) 


(4b) 


where  n =0,1,2...;  k = 1,  2,  3, . . . ; a = 1,2.  With  use  of  the  recurrence 
relations  (4a),  the  series  (3)  may  be  telescoped  as  follows 


(U,of’k|y“’k,.t,  = C' 

+ Jr«*.!Osta.k>[r 

(a»k)  (a,k)  c 

where  C * , S are  formal  differential  operators  defined  by 

Ja,k)  , . ( a)  1 (C£,k) 

C = cosh(c  y a ) , 

Q(o:»k)  (a)  1 (a,k)  -1  . (a)  1 (a,k) 

S = (c  y 3t)  smh(c  y a ) . 


Ja,k), 


,(flt.k). 


(5) 


(6) 


The  functions  <7(0)  a’  (t)>  u(o/a'^(t)  represent  stress  and  displacement  at 
the  midplane  of  the  (<v,k)  layer;  these  serve  as  dependent  variables  in  the 
ensuing  analysis. 

3.2  Transition  to  a Continuum  Theory 

Substitution  of  the  expansions  (5)  into  the  interface  conditions  (2) 
furnishes  the  following  differential  (in  time)  - finite  difference  (in  space) 
equations  with  dependent  variables  defined  at  discrete  points  along  the 
y-axis: 

1 

ST'(cr 


J(a),„(a,k+l)1..(o,k)1  _(P)(f3,k)  (a)cr(a)"  * (a),  (a,k+l)  Ja,k) 


C (u 

d 


(a,  > - 2c  vr  - h'“'E'  - -to) 

Jo),  (a,k+ 1 1 (o.kl  ,_(P)  (P,k)  (OHOUO I 2,  (o,k+l)  (o,k) 

' to(0)  <0(0)  > - 20  °(0)  ■ P 11  S ' ^0)  ) * 0 • 


%)  1 * °- 


(7) 


44 


W 


a,  P = 1. 2;  or  / p. 


If  one  now  defines  the  fields 


(u.c},aV.t)  « (u.af-'V)  (or  y€(*‘°-kLh(«).  Y(“'k|+h,«l)  . 
.(0.0)  (or  iP-W)  . 


a,  P = 1,2;  a ^ P;  and  if  one  notes  that 
(u,c);“;k,(.,  = (u,o)<“,(YI“;klt) 


then  (7)  can  be  rewritten  in  the  form 

Joe)  Joe).  . (a).  , 7JP)(P\  .. 

C [u  (y+A,t)  + u (y-A,t)]-2C  u (y,t) 


Aa)~(ot)  (a)  (or).  ..  (a)  . . n _ n 

h e S [a  (y+A,t)  - a (y-A,t)j  = 0 


JarL  (or),  ,A  ..  , Jar),  ...  jJP)(P),  . 
C X ct  (y+A»t)  + a (y-A,t)]  - 2C  a (y.t) 


(«).(a)Ja)^2  (a)  (a) 

-p  h S stfu  (y+A»t)  - u ( y- A»t)3  =0 


(10a) 


(10b) 


0 i,P  = 1,2;  a / P-  Equations  (10a,  b)  were  obtained  by  replacing  the  fields 
0<P)(y(0,k),  t)  and  Y^P>K)  + A,  t),  etc.,  by  the  smooth  fields  o(P)(y,  t), 

J“)(  v+A,  t),  etc.,  defined  over  the  entire  domain  of  the  composite;  here 
y{or,K)  denotes  the  midplane  of  the  (o,k)  layer  (see  Fig.  1). 

The  functional  difference  equations  (10a,  b)  contain,  as  a subset, 
the  solutions  of  the  finite  difference  equations  (7).  That  is,  the  new 
fields  satisfy  (7)  at  y = Y(P»10.  Therefore,  the  new  dependent  variables 
v/cr),  <j(a)  are  functions  defined  for  all  y which  assume  exact  values  at 
layer  midplanes.  Assuming  the  admissibility  of  such  expansions,  all 
sums  and  differences  are  now  expanded  about  A = 0 in  a Taylor  series; 
this  leads  to  the  partial  differential  equations 

da) c Ja)  - t!a)LEW'lda%  - Cmu,m  = 0 . 

y 

joe)-  (a)  .(or)  (a) Jar)-  2 (or)  J(3)  (P) 

C C<7  - h AP  S Sg^?tu  - C a = 0 , (11) 

or,p  = 1,  2;  or  / p,  where  C,  S are  formal  differential  operators  of  the  form 

C = cosh(A3  ) , S = (AS  ) sinh(A5  ) • (12) 

y y y 

Equations  (11)  constitute  a continuum  model  of  the  composite.  The 
dependent  variables  o(a)  of  (11)  can  be  used  to  reconstruct  the 

microstructural  fields  through  (5)(noting  (8)). 


3.3  Scaling 

Let  the  typical  ' macrodimens  ion"  of  the  problem  be  £ and  the 
typical  "microdimens ion  " be  A.  The  quantity  f represent  a dominant 
signal  wavelength  associated  with  wave  motion.  Introduce  the  non- 
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dimensional  quantities 


. . . Ac r)  _ (or).  2 (a)  _ (or), 

y/i  , t £ tc / f , c e &/*•  , I = a /pc  , T)  = u 1 1 

(13) 

where  p,  c denote  appropriate  mixture  density  and  wave  speed.  Then  (11) 
becomes 


W . o 


where 


Jcr]k  Ja)  ~(a)-  (a)  2 h)  ~((3)  (p)  2 (p)  _ 

C S3.T  - P C8  d_7)  - p 8 ST?  =0  » 


C(a)  = cosh(fy°^3  ) » S(a)  s (Cya\  ) ^inhtcy®^  ) , 

T T T 

C £ cosh(fB^)  , 5 = (CS^)  1 sinh(f3^)  » £ (c/c^ln^ 

(or)  (or)  —(or)  (or)  (or).  (or)  2 (a)  ~(or) 
n t n /i  , p = n p !p  , l £ pc  n /E 


(15a) 


(l  5b) 


(15c) 


The  governing  equations  (14)  are  now  scaled  such  that,  if  y£(0,  t) 
and  T€(0,  l/c),  then  £e(0,  1)  and  Tf  (0,  1);  that  is,  the  macrodimension  is 
now  0(1)  while  the  microdimension  is  O(C).  Thus,  it  is  evident  that  all 
foregoing  Taylor  expansions  about  A=0  are  in  fact  expansions  about  C = 0; 
further,  all  operators  are  power  series  in  e(2 Q> , n =0,1,2.  ..  . Such 

expansions  are  not  required  to  converge,  but  only  be  asymptotic  as 
f -♦  0 [4]. 


3.4  Binary  Mixture  Theories 

With  some  algebraic  manipulations,  equations  (14)  may  be  re- 
written in  a form  characteristic  of  binary  - type  mixture  theories  as 
follows: 


. r.tJ1*).  ~.(or)  (a)-  2 I na+1o 

S^CC^C  - P 8 SoTr)  j = (-1)  ft  » 

. r-Ja).  (»)  >)>):>).  _ ..ar+l 

S^IC-C  - Z S SL  j - (-1)  P » 

a = 1,  2.  Here 


2 _ J 2)  (2) 
f P £ C I 


(1)  (1)  2 (2)  (2)  Cl)  (1) 

CL  ,<ft£Crj  - C rj 


C*  £ 1/21  + (c&c)2/41  + Ua?)4/6l... 


(16a) 


(16b) 


(1 6c,d) 
(I6e) 


The  bracket  in  (16a)  suggests  a momentum  equation  while  that  of 
(16b)  has  the  character  of  a constitutive  relation.  The  quantities  P and  ft 
represent  constituent  interaction  terms. 

A mixture  theory  of  order  n is  now  defined  as  those  equations 
obtained  from  (16a,b)  by  truncating  the  operators  and  as 

well  as  the  in  P and  ft,  after  0(c2n“2)  terms.  In  particular,  a 
mixture  theory  of  order  1 is,  under  this  definition,  given  by 
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aJs„r/2  - pa2r?)<Qf)  = (-i>®+1p*  . a„(a„r)/2-eE)(0f)  = (-i)0+Ia*  d?a,b) 


, V . jf2)-r 


(1)  2 (2)  (I) 

» « R-'  s r?  - r? 


(1 8c,d) 


Equations  (17a)  may  be  interpreted  as  mixture  momentum  relations  and 
(17b)  as  mixture  constitutive  equations.  The  interaction  terms  in  this 
formulation  are  proportional  to  displacement  and  stress  differences, 
respectively. 


The  operator  preceding  the  left  side  of  each  equation  (17a,b) 
can  be  eliminated  by  expanding  the  dependent  variables  in  regular 
asymptotic  series  in  t.  The  result.upon  retaining  0(<  ) terms  only,  is 


(»p>  . (op)  2 <«)  , ,.!,(»)- 
C T 


- 2 ( i ) (2) 

p = -rj  ) 

T 


(18) 


where  £ is  a constant  which  is  a complicated  function  of  composite 
geometry  and  material  properties.  The  quantities  p(aP)  here  play 

the  role  of  "partial"  stresses  and  densities;  their  physical  interpretation, 
however,  is  not  obvious.  Details  concerning  the  derivation  of  (18), 
which  shall  be  called  a modified  (first  order)  mixture  theory,  can  he 


found  in  [41. 


Let  us  now  consider  the  limit  as  c -*  0 in  (16).  If  it  is  assumed, 
that  t and  c may  be  selected  such  that  £°(  ) =0(1),  3»(  ) =0(1)  , 
n = 1,2...,  where  ( ) represents  7 or  as  0(i.  e.  , as  the  ratio 

of  micro-to  macro-dimensions  vanishes),  then  (16)  reduces  in  dimen- 
sional form  to 


2 (1)  (1)  (2)  (2)  (DMl)"1  (2)~(2f2  2..  («)  .... 

ray  - (n  p +n  p )(n  E +n  E )8tJlu*crj  =0.  (19) 

Therefore,  the  theory  is  nondispersive  as  c -»  0,  a result  to  be  expected. 
Equation  (19)  implies  that  the  appropriate  mixture  density  p,  modulus 
E,  and  wave  speed  c,  are 

(!)(!)  ( 2)  (2)  ~-l  (Dcrdf1,  (2)~(2)“ 1 T>~ 

p = n p +n  p , E = n E +n  E , c =E/p.  (20) 


3.5  Phase  Velocity  Spectra 

It  is  evident  that  the  foregoing  theory  represents  a long  wave- 
length, low  frequency  approximation  of  composite  wave  motion.  A 
natural  question  concerns  the  range  of  validity  of  this  approximation.  A 
convenient  indication  of  the  comparative  accuracy  of  an  approximate 
theory  can  be  obtained  by  studying  the  phase  velocity  spectrum. 

Assuming  steady-state  sinusodal  wave  trains,  equations  (16),  (17)  or  (18) 
become  a set  of  homogeneous  algebraic  equations  from  which  the  relation 
between  wave  number  k and  phase  velocity  Cp  is  calculated.  A typical 
comparison  of  exact,  first  order,  and  modified  (first  order)  mixture  phase 
velocity  spectra  is  shown  in  Fig.  2 for  a Polymethyl  methacrylate 
(PMMA)  - 304  stainless  steel  composite.  The  data  shown  represents 
the  first  mode  of  propagation  down  to  the  cutoff  frequency  associated 
with  the  first  "stop  band".  The  exact  spectrum  results  from  retaining 
all  terms  in  the  operators  of  (16)  and  is  identical  to  that  given  by 
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Rytov  [221.  Comparisons  such  as  Fig.  2 indicate  that  reasonable 
accuracy  may  be  expected  from  a theory  of  order  I down  to  wavelengths 
equal  to  several  microdimensions  for  problems  where  the  energy  is 
partitioned  primarily  in  the  first  mode.  This  is  quite  a remarkable 
result  for  such  an  elementary  theory.  Additional  details  concerning  the 
phase  velocity  spectrum  as  well  as  transient  pulse  calculations  and  an 
associated  discussion  concerning  initial  and  boundary  conditions- can  be 
found  in  [4-9]. 

4.  EXAMPLE  2:  HEAT  CONDUCTION 


For  the  purpose  of  illustrating  the  Regular  Asymptotic  Method, 
consider  a linear  thermal  diffusion  process  in  a periodic,  two-dimensional 
array  of  unidirectional  fibers  of  arbitrary  cross  section  embedded  in  a 
matrix,  as  illustrated  in  Fig.  3a.  Let  a "cell"  be  associated  with  each 
fiber  as  depicted  in  Fig.  3b.  Each  such  cell  consists  of  regions  A^  and 
occupied  by  the  fiber  and  matrix,  respectively,  with  unit  outer 
normals*  N^^  and  N^>.  The  interface  between  the  two  constituents 
shall  be  denoted  by  I,  and  the  outer  boundary  of  the  cell  by  C.  With 
respect  to  rectangular  Cartesian  coordinates  xj,  x^,  x^  as  shown  in 
Fig.  3a,  let  the  composite  occupy  the  domain  0 £ £ L,  - ® ® , 

- os  < Xj  < = . Assume  that  the  two  constituents  in  each  cell  are  homo- 
geneous and  isotropic  with  a perfect  interface  (no  thermal  resistance). 
Finally,  let  the  initial  conditions,  and  the  boundary  conditions  on 
Xj  = 0,  L be  such  that  the  temperature  field  is  s imilar  in  each  cell.  In 
view  of  the  lastpremise  it  is  sufficient  to  consider  a typical  cell  with 
lero  heat  flux  normal  to  the  boundary  C.  Consequently  the  basic  equations 
for  the  temperature  fields  T^  and  heat  flux  vectors  a = 1,2,  are 

(Q.  . + dT ja)  = 0 , (Q.  + kT,  .)(Qf)  = 0 on  A(#)  ; (21a,b) 

J.  J 11 


<S(ZI 


0 on  C;  T 


,(D 


T(2)  and  (Qp!)  - Qp2))Np2)  =0  on  I. 


(2  Ic,d) 

Equation  (21a)  is  the  conservations  of  energy;  the  relation  (21b)  is  the 
Fourier  heat  conduction  law;  equation  (21c)  is  a symmetry  condition;  and 
(2  Id)  represents  t.ie  interface  continuity  of  temperature  and  normal  heat 
flux  (note  N3^  = 0).  The  superscript  or  = 1,2  refers  to  the  fiber  and  the 
matrix,  respectively.  The  quantities  d,  k denote  heat  capacity  and 
thermal  conductivity,  respectively.  The  notations(  ),[  = &(  )/dX[  and 
(')  e )/»t  have  been  introduced  for  convenience. 


4. 1 Scaling 

Equations  (21),  together  with  initial  conditions  at  t=0  and  appropri- 
ate boundary  data  on  X3  = 0,  L specify  a well  posed  problem  involving 
three  spatial  variables  Xj  and  time,  t.  The  objective  of  the  subsequent 
analysis  is  to  derive  a continuum  mixture  theory  governing  the  macro- 
scopic diffusion  process  which  involves  only  one  spatial  dimension  (x^) 

‘‘''Latin  and  Greek  subscripts  imply  Cartesian  tensors  and  the  usual  sum- 
mation convention  with  ranges  from  1 to  3 and  1 to  2 respectively. 
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and  yet  reflects,  at  least  approximately,  the  effect  of  conduction  on  the 
microscale.  To  this  end,  let  A and  A be  associated  with  typical 
"macroscopic"  and  "microscopic"  observational  dimensions,  respectively. 
These  lengths  may  be  defined  in  terms  of  characteristic  thermal  diffusion 
times  in  the  longitudinal  and  transverse  directions  according  to 


2 2 

t.  A.nd  >A  /k  . t...sd.  .A  /k 

(A)  (m)  (m)  (A)  (m)  (m) 


(22) 


where  d(m)»  k(m)  denote  mixture  heat  capacity  and  thermal  conductivity 
(to  be  defined  later).  In  addition,  it  will  be  convenient  to  introduce  a 
mixture  heat  flux  Q(m)  based  upon  a reference  temperature  T,  and  a 
parameter  f as  follows: 

°(m)*k(m)T/A  ' • (23> 

The  quantity  t represents  the  ratio  of  micro-to-macro-dimens ions  of  the 
compos  ite. 

With  the  aid  of  the  foregoing  definitions,  the  following  non- 
dimensional  variables  are  introduced 

($3,CV  e <x3'xp)/A  , T = , (Dv  CDfl/  ) s (QvQRf  */Q/mi  . 


(A)  w3'  '“P' 


3*  (m)  ' 


>)  J*)/T,  » ,(a)  .(a)  .(a)., 

7 = T /T  , d = d /d,  , ft  =k  /k  , 

(m)  (m) 


(23a) 

(2.3b) 


whence  (21)  become 


(2.  . + d: r)(ff)  = 0 , («  , (Z?Ja)  + = o on  G(°l  (24a,b) 

J.J  3 p 3 p 

(2)  (2)  Jl)  (2)  . (1)  (2)  (1) 

2 v = 0 on  C ; a = 7 , (2„  -Do  )v„  = 0 on  J . (24c, c5 

P P P P P 

(3)  {&) 

In  (24)  p = 1,2  and  G , C denote  A , I,  C in  nondimens  ional  form; 
the  vectors  are  unit  outer  normals  to  the  boundaries  of  (sa^;  partial 
derivatives  are  now  defined  by  ( ),  = 3(  )/d£^»  ( ) =5(  )/dT  . 

4.2  Mixture  Equations 

Binary  mixture  equations  for  this  problem  can  be  immediately 
obtained  by  averaging  (24a)  over  according  to 


3(aa)(£3,  T)  = G( 


M ,n  Jet). 


j * £3’  1^^2 


(or)  (a)  (a) 

where  3 represents  2 or  3 This  yields,  after  a little  algebraic 
manipulation 

9£>  + rv  _ P,c-':3mvn>ds 

J 


r-lp„(2)  (1). 

= c ;dp  vp  ds  , 


(26a,b) 
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where 


(orp)  (B)  (aa) 
J3  “ J3 


. n<or)  _(ar)  (or)  (a)  (1)  (2) 

= n i , n = G /(G  +G  ) . 


(orp)  (ap)  (a) 

The  quantities  2^  , i , n represent  partial  heat  fluxes,  partial 

heat  capacities,  and  volume  fractions,  respectively.  The  variable  p is 
an  interaction  term  reflecting  heat  transfer  from  the  fiber  to  the  matrix. 
At  this  stage  the  mixture  equations  (26)  are  exact.  Approximations. arise 
when  one  attempts  to  model  the  Fourier  expressions  for  heat  conduction, 
and  the  interaction  term. 


4.3  As^ 


itotic  Exoansions 


A fundamental  premise  is  now  introduced:  the  ratio  of  the 
characteristic  thermal  diffusion  times  in  the  transverse  and  longitudinal 
(fiber  axis)  directions  in  small  compared  to  unity,  i.  e.  , 

= • ,28) 

Equation  (28)  is  appropriate  for  many  composites  used  for  thermal  pro- 
tection. The  premise  (28)  suggests  the  following  regular  asymptotic 
expansion  for  all  dependent  variables,  denoted  by  G^: 


g'  Tjf)  = £ « c£\(«..t) 

1 n=0  (2n)  1 


If  (29)  is  substituted  into  the  governing  equations  (24a,  b)  and  thf  coef- 
ficients of  similar  f-order  are  equated,  one  obtains  a system  of  equations 
for  each  n = 0,  1,  2.  . . . In  what  follows,  a mixture  theory  is  developed 

based  on  the  lowest  order  system. 

4.4  Interaction  Term 

The  lowest  order  system  corresponding  to  (24a,  b)  is 


• (™)  ( *7 ) ( rv ) 

(D.  • +2ir),m  = 0 , (2  + -J/rvx  = 0 » (u,  , = 3",  0 = 0)  . 


(30a,b,c) 


Equations  (30c)  yield 

(or)  (or), 
ff(0)  " °(0)(f3’  r)  * 

whence,  to  lowest  order  accuracy,  equations  (25),  (27)  furnish 


(orp)  _ (orp)^(arp) 

3,  — - o,  , (£ ,»  r)  . 


3 '^3*  ' ' * 

In  order  to  complete  the  mixture  formulation,  the  functional 
dependence  of  the  interaction  term  P on  the  averaged  temperatures  7). 
must  be  determined.  For  this  purpose  it  is  necessary  to  consider 

and  to  satisfy  the  continuity  of  temperature  including  O(c^)  terms. 
To  Degin  this  task,  one  finds  from  (32)  and  (24b)  that 


*P(0) 


(or)  Jar) 

’ '(2),  (3  * 


Next,  with  use  of  (33)  and  (30b,  c),  equation  (30a)  becomes 
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W 7(2rpp 

(a) 


= <p(af)(<3.r)  . 


(34) 


The  functions  <3  can  be  related  to  p,  defined  by  (26b),  by  integration  of 
(34)  over  of0^,  use  of  Gauss'  Theorem,  and  the  symmetry  condition  (24c); 
the  result  is: 


(or)  (a)  a 

a o = (-1)  p . 

Thus,  one  obtains  the  following  relations  for  3 


<«>. 

(2)’ 


(or),  (a)  (a)  . . ..or  Jot) 

n (ft  ^(2)'  p ' p P on  G (or,  P = 1,2)  . 

The  appropriate  boundary  conditions  are,  from  (24c,  d), 

$)'  ^2)  * (^o!  - ^oJ,/eZ  0,1  J • 


(35) 


(36) 


(37) 


Equations  (36),  (37)  constitute  a "micro"  boundary  value  problem 

in  the  £j,  - plane;  the  solution  is  unique  within  a function  Tr(£3,  r).Once 

are  known,  one  can  write 
Iw 


where 


3J°f)  = [S&!  + f23C  = ‘?(“J(C3,t)]  + e2p7-(a)(^1,42)  + 0(f4) 


_(D  _ JD  „*(2)  _ J 2)  . ,~(2)  ~d)  , 2 

PT"  “ °(2)  ’ _ °(2)  + °(0)  " 0 (0 ) ^ e * 


Upon  averaging  (39),  one  obtains 


Jaa).  . ~(U  w 2 (aa) 

S (£3*  t)  ~ ?(o/^3’t^  + c 


Hence,  the  interaction  term  can  be  written 

„ _ J2a).,  2 _ (la)  ,(2a).-l 

P = a(7  -o  )/(  , a s (5J*  -5-  ) 


(38) 


(39) 


(40) 


(41a,  b) 


Equations  (4  la),  (32),  (26a)  close  the  mixture  theory.  The  primary 
result  is  expression  (41b),  which  can  be  used  to  determine  the  interaction 
coefficient,  a,  for  arbitrary  fiber  and/or  cell  geometry  via  the  solution 
of  a time  independent  problem  defined  over  the  unit  cell.  In  terms  of 
Z*^aK  this  problem  takes  the  form: 


n(“W«“i.p . i-i)c 


on  C 


(a) 


(42) 


.(2)  (2)  „._,(r)  -i. ( 2 ) Jl)  JD  (2)  .(21 

JT^'p  Vp  =0  on  c ; =Z*'  . (ft  rr-"p  -ft  T*p  JVp  = 0 on  J,  (43) 


7-^=0  at  point  0~C^^ 


(44) 


In  general  it  is  difficult  to  obtain  an  analytical  solution  to  the 
above  micro  boundary  value  problem.  As  an  alternative,  a finite 
element  procedure  has  been  proposed  [21]  which  may  be  utilized  for 
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arbitrary  two-dimensional  cell  geometry.  This  procedure  is  based  upon 
a modified  Reissner-type  variational  principle,  the  corresponding 
functional  of  which  is  defined  by 


n=  I ; 

®- 1 (a)  ,P  ,P 

n 


(-i)' 


.(or) 


(a) 


(45) 


which  is  defined  in  terms  of  J*0  on  G a and  2*  on  J.  The  equations  for 
follow  from  the  requirement  n be  stationary  with  respect  to 
arbitrary  variations  of  and  2*;  here  2*  is  a Lagrange  multiplier 

which  physically  represents  the  heat  flux  normal  to  the  boundary  j. 


4.5  Numerical  Results 


The  finite  element  procedure  has  been  used  for  the  solution  of  the 
temperature  microstructure  problem  and  the  interaction  coefficient  for  a 
variety  of  geometries  and  combination  of  material  properties.  Details  of 
such  calculations  can  be  found  in  [21  j.  Here  several  results  of  particular 
interest  are  noted;  properties  used  are  listed  in  Table  3. 


A frequent  fibrous  composite  geometry  is  a hexagonal  array  of 
circular  cylindrical  fibers.  In  practice  it  is  common  to  approximate  the 
cell  geometry  of  this  case  by  concentric  circular  cylinders  [20], 
Numerical  results,  based  on  the  theory  presented  herein,  related  to  the 
accuracy  of  this  approximation  are  shown  in  Fig.  4,  where  the  inter- 
action coefficient  has  been  calculated  as  a function  of  fiber  volume 
fraction  and  thermal  conductivity.  From  the  data,  the  conclusion  can  be 
drawn  that  for  a practical  range  of  fiber  volume  fractions, the  concentric 
cylinder  approximation  can  be  used  without  any  significant  loss  of  ac- 
curacy. 


For  calculation  of  the  interaction  coefficient,  the  concentric 
circular  cylinders  approximation  based  on  equal  fiber  volume  fraction 
can  also  be  used  for  composites  containing  square  fibers  arranged  in  a 
square  array.  This  is  borne  out  by  the  results  shown  in  Fig.  5. 

The  fields  £3)  nee<3ed  to  calculate  the  interaction  term  can 

also  be  used  for  the  calculation  of  the  temperature  distribution  in  a unit 
cell.  From  (38)-(40)  one  has 


+ at • (46) 

In  the  right  side  of  (46)  the  only  quantity  dependent  upon  inplane  coor- 
dinates and  £2  It  follows,  therefore,  that  curves  of  equal 

are  also  isothermal  lines  within  the  framework  of  the  mixture 
theory.  For  this  reason,  and  to  illustrate  the  type  of  temperature  micro- 
structure  that  can  be  obtained  from  the  mixture  theory,  contours  of  equal 
"^®),  suitably  normalized,  are  given  in  Fig.  6 for  square  fibers  in  a 
similar  array. 

Finally,  an  indication  of  the  accuracy  of  the  proposed  mixture 
model  and  the  construction  procedure  can  be  obtained  by  comparison  of 
mixture-predicted  response  with  2D  or  3D  computer  code  calculations. 
This  has  been  accomplished  in  [20]  for  the  case  of  concentric  circular 
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cylinders  where  two  spatial  dimensions  and  time  are  involved. 

Calculations  were  conducted  to  determine  the  evolution  of  the  temperature 
field  in  a quescent  half-space  £ 0 subject  to  the  boundary  condition 


T(0,  r,  t) 


1,  Ostst 

0 


0,  t > t 

o 


(47) 


Both  the  mixture  equations  and  the  full  2D  problem  were  solved  by  finite 
difference  schemes  [20];  the  latter  was  used  as  the  correlating  norm  for 
estimating  the  accuracy  of  the  mixture  theory  solutions.  Fig.  7 illustrate 
typical  average  temperature  profiles  in  the  two  constituents  after  a short 
time  following  the  termination  of  the  temperature  pulse.  The  ability  of 
the  mixture  theory  to  predict  the  temperature  microstructure  is 
illustrated  in  Fig.  8,  the  radial  distribution  of  temperature  obtained  from 
the  mixture  theory  is  almost  identical  to  the  exact  solution. 


4.6  Remarks 

The  significant  feature  of  the  above  mixture  theory,  and  of  similar 
theories  proposed  for  wave  propagation  [16-19],  is  that  it  converts  what 
is  essentially  a three-dimensional  problem  to  a problem  involving  a 
single  spatial  variable,  without  losing  the  essential  details  of  local  field 
distributions.  The  reduction  in  the  number  of  dependent  variables  leads, 
of  cou.se,  to  a substantial  increase  in  numerical  efficiency. 


5.  SELECTED  APPLICATIONS 

Two  applications  involving  nonlinear  wave  propagation  are  briefly 
reviewed  in  this  section  in  an  effort  to  demonstrate  the  utility  of  mixture 
theories  with  microstructure.  The  first  involves  delamination,  of 
laminated  composites.  The  second  concerns  finite  amplitude  elastic- 
plastic  wave  propagation  in  fibrous  composites. 


5.1  Wave  Propagation  and  Debonding  in  Laminated  Composites 

A "first  order"  mixture  theory  with  microstructure  has  been 
developed  in  [ 12]  for  longitudinal  wave  propagation  and  debonding  in  a 
laminated,  binary  composite  with  periodic  microstructure.  The  case 
considered  concerns  small  isothermal  deformations,  linear  elastic  con- 
stituents, and  a Mohr-Coulomb  interface  failure  and  slip  model.  In  view 
of  the  last  item  the  problem  is  nonlinear.  With  respect  to  Fig.  1,  a 
condition  of  plane  strain  was  assumed  in  the  X3-direction  and  wave 
motion  yielding  symmetric  velocity  Vj  and  anti-symmetric  velocity  y, 
distributions  with  respect  to  the  X2~coordinate  within  each  lamina. 
Theoretical  construction  was  based  on  the  Differential-difference  Method. 
In  dimensional  form,  and  with  reference  to  Fig.  1,  the  relevant  equations 
are 
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0'«p>i.p<«p>s,  = ,.if+1P 


1 


(op)  £ (ay)  (ya) 
°11  = r,c  “1.1 
y=1 


(48a, b) 


P = K(Uj  -u{ 


(la)  • (2a)^ 


if  P<lP 


cr 


A-BO 


p = e 


h(1)+h(2) 


22 up.|Peri 


i 


i 


(48c) 


Equations  (48a)  are  mixture  momentum  equations  while  (48b)  represent 
mixture  constitutive  relations.  The  quantity  P denotes  a momentum 
interaction  term.  The  variables  p(^P^  denote  partial  quantities  and 

are  defined  by  ( = n*a\  where  ( is  an  averaged  variable 

and  n^  = h^/h*1)  + h(2).  The  quantity  |Pcr|  is  related  to  the  interface 


shear  stress  crj^  an^  normal  stress  o by 


■(2’)| 


I = A-Bc*  (a  _ s 0)  » 


Pcrl  " lCT12  ' ~ * 22  ,w22 

cr 


ah  - (n(>y'y2,u.,+n(2y2yl,u(2^/(n(.y2)+n(2)E(i))  _ (49| 


The  constants  A,  B represent  interface  cohesion  and  frictional  coefficients, 
respectively.  The  constants  <?a^,  K are  given  in  terms  of  the  composite 
properties  and  geometry  by 


W .«)/ Ji).,d)u  (2)  <2\  <aa>  - >u«)  _ £ 


(o)2 


= .3fi  /jl  / (h  *rh  ) (n  +n  fi  ) , c = n E 


<*> 

c * E • L “ (1)  (2) 

n n 


, Eia)  = (>+2p)(a),  (o,  (3  = 1,2;  or#), 

(BO) 


where  X>  P are  Lame  ' constants. 

Considerable  insight  into  the  modeling  capability  of  the  above 
mixture  theory  can  be  gained  by  comparison  of  theoretical  calculations 
with  the  experimental  impact  data  on  delaminated  plates  reported  in  [23], 
The  composites  in  this  study  consisted  of  alternating  layers  of  Poly- 
methyl methacrylate  (PMMA,  Rohm  and  Haas  type  A)  0.762  mm  thick, 
and  6061-T6  aluminum  0.792  mm  thick.  The  laminae  of  the  composite 
were  oriented  perpendicular  to  the  impact  plane,  and  struck  by  a 
projectile  from  a 10cm  bore  light  gas  gun.  An  aluminum  buffer  plate 
1.0cm  thick  was  placed  at  the  rear  of  the  composite  to  improve 
planarity  of  the  transmitted  wave  front.  A transparent  Dynasil  100 
window  followed  the  buffer,  and  a thin  mirror  was  vapor  deposited  at  the 
buffer  window  interface.  The  motion  of  the  mirror  was  monitored  to 
within  ±1.5  "*  10*  5mm  by  means  of  a displacement  interferometer.  The 
experimental  configuration  reproduced  from  [23]  is  depicted  in  Fig.  9 . 
The  shot  matrix  is  summarized  in  Table  1. 

The  mixture  theory,  consisting  of  (48)  was  coded  in  finite  dif- 
ference form  and  solved  numerically.  Material  properties  for  the 
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composite  laminates,  flyer  plate  material,  buffer  and  window  materials 
are  listed  in  Table  1,  Input  to  the  computer  code  required  only  the 
material  properties  of  all  constituents,  goemetries  of  the  test  configura- 
tion, and  flyer  impact  velocity.  Numerical  values  for  A and  B were 
determined  by  a parametric  study  on  the  results  of  experimental  shot 
number  2.  They  were  found  to  be  A = 0.01  x 10^  dynes/cm^  and  B = 0.50. 
All  other  correlations  were  performed  using  these  values.  Typical 
results  are  depicted  in  Figs.  10-12.  As  can  be  seen,  agreement  is 
generally  excellent.  Both  theoretical  and  experimental  results  clearly 
exhibit  the  effects  of  bond  breakage  and  delamination.  This  is  manifested 
by  the  appearance  of  a precursor  at  the  front  of  the  wave  profile  and  by 
greatly  reduced  oscillations  behind  the  wave  front  when  compared  with 
both  theoretical  and  experimental  results  for  perfectly  bonded  wave- 
guides (e.  g.  , see  [12]). 


5.2  Nonlinear  Wave  Propagation  in  Fibrous  Composites 

A binary,  first-order  mixture  theory  with  microstructure  has 
been  developed  in  [ 17],  using  the  Regular  Asymptotic  Method  and  the 
concentric  cylinders  approximation,  for  nonlinear  waveguide-type 
propagation  parallel  to  the  fibers  of  a unidirectionally  reinforced  fibrous 
composite  with  initial  periodic  hexagonal  array.  Fig.  13.  The  resulting 
model  incorporates  the  effects  of  thermodynamics,  finite  deformation, 
and- elastic-plastic  constituents.  The  constitutive  relations  for  this 
analysis  consist  of  a Mie- Griineise"  equation  of  state  which  relates  mean 
pressure,  density,  and  internal  energy,  and  a von  Mises-type  yield 
^cirteria  and  associated  (J2)  flow  rule  which  govern  the  stress  derivators. 
Under  the  assumption  of  an  adiabatic  process  and  isotropic,  homogeneous 
constituents,  the  following  mixture  theory  was  obtained  for  composites 
with  perfect  bonds 

(a)  Continuity: 

(2)  (3 a)  /ora)  (ora) 

L n i =0  , f s — T(p  v 

a = l 


(ora)  , • (ora)  . (ora) 
3 ).3*P  ]/ P 


(b)  Momentum: 
;33 


P (Of),  ,c  >(aa)  (a)  (ora)  (oa)  (ora)  nar+l 

[n  (-p+S,,)  ],  3-n  p Dt  v3  =(-l)  P , 


where 

| nMB{aVai) 

0r=l 


U,,  , «o,  ,c  ,(“a)  (or)' 
P+(-l)  (-p+S33)  n, 3 


(or) 


o,  Ua)  (2a). 

, ,.cr  (aUflliaa)  8<V3  "V3  } 

- (-1)  n B f,  


(51a) 


(51b) 


I,M)(  ),  ),  D<°a)<  ),(Vv^a)(  ).3,  (51c) 
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(c)  Energy: 


(2/n(3)2|(i.2!1«)t|„(.)2+lnn<'))  . 


Of  = 1 

_ , ,c  >(Ia)  (D.da) 

R = (-p+S  ) n £ 

rr 

(d)  Caloric  equation  of  state  (Mie- Grvineisen); 

(oa)  (aa)  (aa)  (aa)  (aa) 

p = g(p  )+r(p  )p  e 

(e)  Stress  derivator  constitutive  relations: 

(oa)_(aa)  4 (aa)  lf(aa)  Jla).  -(la)  _ J-(la) 

X S33  =l(V3,3-2f  ’ S99  " S~  "“’S”  ’ 


rr 


(2a)  (2a)  2.  (aa)  1 (aa)  3 1 

* S-  '•?',3.3-I£  (‘  + (2)ln  (I)’1  ’ 

n n 


rr 


s'2;*1.  - (s'2a,+s'2*,»  , 

S0  rr  33 


where 


0 if  or  D<“a)j!,“a|  < 0 . 

* 2 3 t 2 

= 3 Y(a)2rs(oa)v{aa)+s(aa)f(aa)  .£  D(oa)j(aa)  ^ Q 
2 33  3,  3 rr  J t 2 


and  J 


(aa)  Y 


.(a) 


2 3 

(f)  Interface  stress  match: 

(-p+S  )(la)(x,,t)  = (-p+S  / a'(x  ,t)  . 

rr  3 rr  3 

(g)  Volume  fractions: 

(la)  (1)  (1)  (la)  (1),  (2)  _ 

n =n£  , n +n  -l. 

(h)  Definition  of  averages: 

JD 


( y 
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irr 
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2 0 


r (1) 

2Trr(  Y dr  , 


(51  d) 

(52a) 

(52b) 

(53) 

(54a) 

(54b) 

(54c) 


(54d) 


(55) 


(56) 


(57a) 
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(2a) 


(2) 


ir(r 


(2) 


-r(1)  ) r(1) 


The  dependent  variables  above  number  14.  These  are  n , p # 
p(aa-\  S^a)  and  e(aa)(o  = l , 2).  The  foregoing  represent  volume 
fraction, density,  velocity  component  in  the  X3~direction,  stress  deviator 
component  in  the  X3~direction,  and  internal  energy,  respectively. 
Equations  (51)  — (56)  constitute  14  independent  relations  for  the  above 
unknowns.  The  independent  variables,  which  were  originally  r,  9,xj,  t 
(see  Fig.  13),  are  now  x^.t  only.  All  equations  are  written  in  an 
Eulerian  description.  It  should  be  noted  that  a = 1 refers  to  the  fiber 
while  a -2  denotes  the  matrix. 


In  an  effort  to  test  the  validity  of  the  first-order  mixture  theory,  a 
comparison  was  made  between  calculations  performed  using  the  quasi- 
one- dimensional  mixture  equations,  and  a formal  two-dimensional 
solution  of  the  concentric  cylinders  problem.  The  mixture  solutions  were 
obtained  by  writing  (50-56)  in  finite  difference  form  and  solving 
simultaneously  for  the  14  dependent  variables.  The  necessary  numerical 
analysis  and  the  resulting  code  is  described  in  [24],  The  two- 
dimensional  solutions  were  generated  using  a well-known  Lagrangian 
finite  difference  computer  code  called  CRAM  [25].  A typical  cell  of  a 
quartz-phenolic  composite  occupying  the  half- space  X3  < 0 was  selected 
for  study.  For  computation  in  the  range  0 s p s 30  K bar  a Mie- 
Griineisen  relation  of  the  form 


(aa)  3 (a)  (a)  (or)  (aa)  (aa) 

p = r H.  TJ  +G  p e 

i=l 


(or)  (ora),  (a)  . 

n * p 4>(I)- 1 


(58) 


was  used,  where  denotes  the  initial  density  of  material  or,  and  , 
G^“)  are  material  constants.  The  constituent  material  properties 
employed  in  all  calculations  are  given  in  Table  2.  Quiescent  initial 
conditions  were  assumed.  Boundary  conditions  consisted  of  a particle 
velocity  of  0.5cm/p.  sec  of  1.0-p.  sec  duration  applied  uniformly  to  both 
constituents  at  x^  = 0. 


A comparison  of  mixture  theory  and  two-dimensional  finite  dif- 
ference code  predictions  of  averaged  constituent  particle  velocity  as  a 
function  of  time  is  illustrated  in  Figs.  14,  15.  The  two-dimensional  code 
prediction,  which  resolves  the  particle  velocity  distributions  in  the 
radial  direction,  where  averaged  according  to  (57)  in  order  to  provide 
direct  comparisons  with  the  mixture  theory.  Results  are  shown  for  a 
propagation  (x3~direction)  distance  of  0.5cm.  Agreement  between  the 
finite  difference  and  mixture  theory  calculations  is  judged  to  be  excellent. 
Finally,  it  is  noted  that  on  the  same  computer  (aUNIVAC  1108),  the  CRAM 
calculation  required  roughly  500  times  the  computer  time  needed  for  the 
equivalent  mixture  calculation. 
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Table  1.  The  Experimental  Configurations. 


Composite 

thickness 

(cm) 

Projectile 

Experiment 

No. 

Material 

Velocity 
(cm/  psec) 

Thicknes  s 
(cm) 

1 + 

0.  820 

Aluminum 

0.001355 

1.  634 

2 

0.776 

Aluminum 

0. 001420 

1.  631 

3 

0.  806 

Aluminum 

0.001144 

1.  571 

4 

0.  803 

PMMA 

0.0021 18 

0.  691 

5 

0.  805 

PMMA 

0.003078 

0.  694 

6 

0.  812 

Tungsten 

Carbide 

0.001289 

0.  975 

7 

0.  810 

Aluminum 

0.001330 

0.656 

8 

0.  809 

Aluminum 

0.001066 

0.  164 

* In  all  experiments  the  buffer  plate  was  1 cm  thick.  The  bond 
thicknesses  were  all  less  than  0.0002  cm  thick. 

^ Experiment  1 exhibited  excessive  flyer  tilt  and  Drumheller  and 
Lundergan  considered  the  data  useless. 


61 


Table  2 


r 

Material  Properties  used  in  TINC  and  CRAM  Calculations 


Material  1 Material  2 


Initial  Density  (g/cc) 

1.  85 

1.29 

Shear  Modulus,  p (dynes/cm2) 

106.0  (109) 

30. 3 (10s) 

Bulk  Modulus,  (dynes/cm2) 

397.6  (109) 

96.0  (109) 

dynes /cm2) 

38.0  (1011) 

0 

(dynes/cm2) 

3 

275.0  (1011) 

11. 3 (1011) 

Mie-C-riineisen  Ratio,  G^°^ 

. 32 

. 52 

Yield  Stress  in  Simple  Tension 

(dynes/ cm2) 

2.  7 5 ( 109  ) 

1.  33  (109) 

Radial  Fiber  Radius,  r (cm)  0.08 

Outer  Cylinder  Radius  r (cm)  0.  20 

a 
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Fig.  10 


